cond-mat/0007395 



DFUB/00-12 
September, 2000 



Bose-Einstein condensation in the 
presence of an impurity 



Paola Giacconi", Fabio Maltoni , and Roberto Soldati" 

° Dipartimento di Fisica "A. Righi", Universita di Bologna, and 
Istituto Nazionale di Fisica Nucleare, Sezione di Bologna, Italia 

^ Department of Physics, University of Illinois at Urbana-Champaign, 

Urbana, IL 61801 USA 



Abstract 

It is shown that Bose-Einstein condensation occurs for an ideal gas in two spatial dimensions in 
the presence of one impurity which is described quantum mechanically in terms of a point-like 
vortex and a contact interaction. This model is exactly solvable and embodies as a special case 
the analogous problem in three spatial dimensions. 
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1 Introduction 



It is well known that Bose-Einstein condensation, a first order phase transition in momentum 
space, can not occur for an ideal gas of free particles in two dimensions. In this short note we shall 
show that, by contrast, the introduction of one point-like impurity on the plane just allows the Bose- 
Einstein condensation to take place. Furthermore, the general pattern of the above phenomenon 
will be obtained from a generating Hamiltonian, which encodes a two parameters model to describe 
the point-like impurity in two and three spatial dimensions. The key point to be gathered is how 
to treat in quantum mechanics the presence of a point-like, or J-like, or null support impurity. 
Formal manipulations involving some kinds of regularized 5-like potentials might drive, in general, 
to misleading and incorrect conclusions, as the latter ones are mathematically ill-defined. The correct 
quantum mechanical framework to treat point-like impurities Q is by means of the analysis of the 
self-adjoint extensions of the Hamiltonian operator which, in the present case, turns out to be just 
a symmetric operator. Consequently, its domain has to be suitably defined in order to obtain a 
self-adjoint Hamiltonian operator which admits a complete orthonormal set of eigenstates. This 
construction will be referred to in the sequel as the inclusion of contact interaction. To be definite, 
let us consider as a starting point of our analysis the following classical one-particle Hamiltonian in 
two spatial dimensions: namely. 



H((j)) = — 



P - -A(r) 

c 



p, r e r2 



(1.1) 



Aj{xi,X2) = —ejk— , r = Jxf + x^ 



2tt '"r^ 

Here the Aharonov-Bohm type vector potential corresponds to the presence of a (5-like vortex 
of flux (/>, which provides a good classical description of one point-like impurity, as we shall better 
specify below. Quantization of the classical Hamiltonian (1.1) leads to a symmetric operator and, 
consequently, one has to face the problem of finding all its self-adjoint extensions. The most general 
solution has been recently obtained in Ref. [Q] and it consists in a four parameter family. However, 
to our aim, we can restrict ourselves to the one parameter sub-family of the 0(2) rotational invariant 
self-adjoint Hamiltonians ||5|, |^. The corresponding spectral decompositions read 

H{a,E^)=Y, / dk ^\l,k){k,l\+'&{-Eo)\il^B){i^B\ , 

l=—oo ^ 

(1.2) 

a = ^ G] - 1, 0] , -oo < Eq < +00 
nc 

t9 being the usual Heaviside's step distribution, in terms of the eigenfunctions 

(r, d\l, k) = ^^Pg.^} ^,(fe, r;a,Eo), k>0, (1.3) 



where the improper eigenfunctions belonging to the continuous part of the spectrum are given by 

^l{k,r) = VkJ\i+^\{kr) , Z e Z - {0} , (1.4) 
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Mk,r;Eo) = A{k;a,Eo)J\^\{kr) + B{k;a,Eo)N\^\{kr) , (1.5) 

in which 

B{k;a,EQ) sin(7ra) 



A{k;a,Eo) cos{7ra) + sgn{Eo) {h'^k^ /2m\Eo\)" ' 
whereas the normaUzable bound state is provided by 

, , , . , , , K /sin(7ra) , . . / tttt 

{rlipB) ='ipB{K,r) = -\ Ka[Kr) , Hk = y/2m\Eo\ . 

IT \ a 



(1.6) 



(1.7) 



Some remarks are now in order. First, as previously noticed, the above spectral decompositions 
precisely provides the correct mathematical framework to introduce and properly describe contact 
interaction in quantum mechanics. As a matter of fact, it turns out that the rotational invariant 
self-adjoint Hamiltonian operators H{a,Eo) do represent a one parameter family, which is labeled 
by the energy scale Eq. In the range — oo < Eq < 0, and only within this range of values, a bound 
state \iPb) exists, whose energy is just Eq. More generally, the physical meaning of the characteristic 
energy scale Eq is given by the resonance energy, according to the following pattern: namely, 

_ / |^o|(sec7rQ)i/l"l ^ if < |a| < 1/2 , < ; , . 

" I \Eo\, in/2 < |a| < 1 , ii;o < ; ^^'^^ 

E,es = £;o|sec7ra|^/l'^l , if 1/2 < |a| < 1 , > . (1.9) 

A further observation is that only the non-integer part of the vortex flux parameter a is actually 
observable, as its integer part can always be gauged away by a single valued phase transformation. To 
sum up, we can say that a general correct quantum mechanical description of one point-like impurity 
is provided by the two parameters family of self-adjoint Hamiltonians of Eq. (^]^). The existence 
of the contact interaction just corresponds to the presence of a specific locally square integrable 



singularity of the wave function at the impurity position - see Eq. (1.5). In the limit Eq —oo 
contact interaction is removed, the domain of the Hamiltonian is that of the regular wave functions 
on the whole plane and the impurity is described in terms of a pure Aharonov-Bohm vortex of 
non- integer vorticity a. If we further take the limit a | 0, i.e., also the Aharonov-Bohm interaction 
is turned off, the two dimensional free particle Hamiltonian is truly recovered (Priedrichs' limit). 
As we shall discuss in the sequel, it is curious that just in the Friedrichs' limit the Bose-Einstein 
condensation disappears in the two spatial dimensional case because, in the presence of contact 
interaction and no matter how weak it is, a non-vanishing critical temperature for the Bose-Einstein 
transition always exists. 

2 One-particle partition function 

In order to discuss Bose-Einstein condensation, it is necessary to compute the average number of 
particles at thermal equilibrium. To this aim, let us first evaluate the diagonal Heat-Kernel and 



the one-particle partition function. According to the spectral decomposition of Eq. (1.2) and after 
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separation of the truly free particle Hamiltonian Hq = H{0,—oo) contribution, it is not difficult to 
verify that the diagonal Heat-Kernel can be cast in the following form: namely, 

G{a,p,Eo;r) = Gint(a, A ^o; + Go(/?) 

= {r\[exp{-f3H{a,Eo)} - exp{-PHo}]\ r) + A^^ 

= /(a; r) + I{-a; r) - 2/(0; r) - /o(a; r) - /o(-q; r) + Io(0; r) 

+ i9(-So)e-^^"|V'B('^r)|2+Xo(a,^o;0 + Ar' , (2-1) 

where the translation invariant free particle diagonal Heat-Kernel is nothing but the inverse square 
thermal wavelength Xt = {h/^/2-KmkT) and we have set 

'^ = ['^ ke-^^'^'l^- f; {J,^^ikr)\^ , (2.2) 



i=0 



Io{a;r) = / — fce"^^ Z^™ [Ukrf , (2.3) 

2( E-r^-T— fcexp{-/3;t^fcV2m} 

" ' °' Jq 27r 1 -|- tan^[7r^(A;)] -|- 2 tan[7r^(A;)] cos (avr) 

X {tan2[7r/i(A;)] J\{kr) + Jlikr) + 2ian[T:^i{k)] J-a{kr) J„(/cr)} , (2.4) 
"2m|/;o| 



tan[TTn{k)] = sgn(E'( 



(2.5) 



h^k-^ 

Now, it is very important to realize that the impurity interaction part of the diagonal Heat- 
Kernel Giotiot, f3, Eo;r) is integrable on the whole plane. This leads to the following result for the 
one-particle partition function 
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asin(7ra) ("^ dx sgn(/;o) e-^l^ol^ 

^ TT Jo 1 + 2sgn{Eo)x\<^\ cos(7ra) + x^l^l ' 

where, as usual, we have denoted by A the area divergence, due to the presence of the translation 
invariant part of the free one-particle Heat-Kernel. The above expression for the one-particle par- 
tition function can be used as a generating form which encodes different specific notable cases. In 
particular, the one-particle partition function in two spatial dimensions and in the presence of pure 
contact interaction can be obtained in the limit a | and reads 

A_ f^dE e-f"^ 

A 



where 



+ u{P\Eb\) , Eb<0, {2.i 

Ar, 
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Notice that in two spatial dimensions the bound state is always present for any —oo < Eb < 0. 



Another distinguished case that can be read off the basic formula (2^) is the three dimensions 
one-particle partition function in the presence of contact interaction. As a matter of fact, thanks 
to dimensional transmutation ||5|, the latter case just corresponds to the value a = —1/2, up to a 
suitable redefinition of the free part: namely, 

V . 1 



Z3dW, Eo) = j^ + ^(-i?o)e^l''°l + -sgn(i^o)e''l''"lerfc(7;p^^) . (2.10) 



3 Results and Discussion 

Now we are ready to discuss the Bose-Einstein condensation for an ideal gas of particles in the pres- 
ence of one point-like impurity, as generally described by the one-particle Hamiltonian of Eq. ( |1.2| ). 
According to the general form ( |2.6| ) of the one-particle partition function, it turns out that the av- 
erage particles density at thermal equilibrium in two spatial dimensions and in the presence of one 
point-like impurity is given by 



2A(l-z) ' "'A{zQ-z) 

where we have set zq = expj/JE'o} and gi{z) = — ln(l — z), whereas 

f TP 1 1? h a sin(7ra)ii^'"l~-^ 

" vrl^ol" + l^oPl"! + 2sgn(^o)(^|^o|)l"l cos(7ra)] ' ^ ' 

It is important to realize that if < the range of the fugacity is < z < zq < 1, whilst < z < 1 
if Eq > 0. Moreover, it is not difficult to prove that, thanks to analytic continuation, the very last 



term in Eq. (3J) admits a finite limit when z ] 1 and Eq > 0. From the above expression (|3.l| ) for 
the average particle density in two spatial dimensions, it appears to be manifest that Bose-Einstein 
condensation occurs only in the presence of the bound state, i.e., only for the sub-family of the 
self-adjoint extensions of the symmetric Hamiltonian ( |l.lD in the range — oo < £"0 < 0. In those 
cases, the critical temperature and/or specific area can be obtained as the unique solutions of the 
equation 



In l-e^-« =-:^(nU . (3.3) 



h^f3 
2TTm 

The three spatial dimensional case can be handled in a quite similar way, as it essentially corresponds 
to the specific value a = —1/2 in Eq. (|3.1j), up to terms irrelevant in the thermodynamic limit: 
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namely, 



+ sgn(^o)i: r''E '''^:" = -'^f±f\ (3.4, 



V Jq 1- zexp{-PE} 

The above equation clearly indicates that Bose-Einstein condensation always takes place in three 
spatial dimensions, in the presence as well as in the absence of the impurity. Nonetheless, the 
actual values of the critical temperature and/or density do depend upon the sign of the parameter 
characterizing the self-adjoint extension of the Hamiltonian. In fact, for Eq > 0, i.e. in the absence 
of the bound state, the critical values are the usual ones as given by the solution of the equation 
Xj, {n)^Y) — C(3/2). At variance, when — oo < Eq < 0, i.e. in the presence of the bound state, the 
critical values can be read off the equation 

gsizo) = A|^(n)3j3 . (3.5) 

The case of pure contact interaction in two spatial dimensions can also be obtained from the basic 
formula ( |3.1| ) taking the limit a ] and treating separately the cases Eq > and — oo < Eq < 0. As 
a matter of fact, the result is 

dE e-^^ 1 

~¥ 1 - 2exp{-/3^} \n^{-E/EQ) + t:^ 

which shows that condensation does not occur when Eq > 0, whereas it appears if — oo < i5^o < 0, 
the critical temperature and/or specific volume being always determined by Eq. ( |3.3D which does 
not depend upon a. In this latter case, the very same formula can be also obtained directly from 
Eq. (pl8|) as it does. 

In conclusion, we have shown in this note that the presence of contact interaction makes it 
possible the occurrence of Bose-Einstein condensation in two spatial dimensions. This phenomenon 
is connected to the presence of a bound state in the spectrum of the self-adjoint Hamiltonian. It turns 
out to be remarkable that the latter circumstance is always there in the pure contact interaction 
case, that means without Aharonov-Bohm vortex interaction. It is in fact worthwhile to notice 
that, in the presence of a non- vanishing vorticity a, the half- family without bound state of the 
self-adjoint extensions does not allow Bose-Einstein condensation, whilst the remaining half-family 
with bound state leads to Bose-Einstein condensation though the critical temperature is independent 
from a. Accordingly, only when — oo < £"0 < a non-vanishing and vorticity-independent critical 
temperature is allowed in the two spatial dimensional case - see Eq. ( |3.3| )- whereas the critical 
temperature deviates from its conventional value in the three spatial dimensional case - see Eq. (13.51) . 



References 



5 



[1] K. Huang, Statistical Mechanics, John Wiley & Sons, Singapore (1987); 
R. K. Pathria, Statistical Mechanics, Pergamon Press, Oxford (1972). 

[2] S. Albeverio, F. Gesztesy, R. Hoegh-Krohn and H. Holden, Solvable Models in Quantum Me- 
chanics, Springer- Verlag, New York (1988). 

[3] Y. Aharonov and D. Bohm, Phys. Rev. 115 (1959) 485-491. 

[4] R. Adami and A. Teta, Lett. Math. Phys. 43 (1998) 43-53. 

[5] C. Manuel and R. Tarrach, Phys. Lett. B268 (1991) 222-226; 

[6] P. Giacconi, F. Maltoni and R. Soldati, Phys. Lett. B441 (1998) 257-264; 
Phys. Rev. D53 (1996) 952-959. 

[7] P. Giacconi, F. Maltoni and R. Soldati, Phys. Rev. B53 (1996) 10065-10071. 

[8] A. P. Prudnikov, Yu. A. Brychkov and O. L Marichev, Integrals and series, Gordon and Breach 
Science Publishers, New York (1992). 



6 



